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^ : ABSTRACT 

■ We present r-^2"-A:" aperture synthesis observations of the circumstellar disk 

surrounding the nearby young star TW Hya in the CO J=2-l and J=3-2 lines and 
associated dust continuum obtained with the partially completed SubmiUimeter 
Array. The extent and peak flux of the 230 and 345 CHz dust emission follow 
closely the predictions of the irradiated accretion disk model of Calvet et al. 
^ ■ (2002). The resolved molecular line emission extends to a radius of at least 200 

^ ■ AU, the full extent of the disk visible in scattered light, and shows a clear pattern 

O ■ of Keplerian rotation. Comparison of the images with 2D Monte Carlo models 



> 



O 



constrains the disk inclination angle to 7°±1°. The CO emission is optically thick 
in both lines, and the kinetic temperature in the line formation region is ^20 K. 
Substantial CO depletion, by an order of magnitude or more from canonical dark 



. cloud values, is required to explain the characteristics of the line emission. 
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Subject headings: stars: individual (TW Hya) — stars: circumstellar matter — 
planetary systems: protoplanetary disks — radio lines: stars — ISM: molecules 

1. Introduction 

TW Hya is the closest known classical T Tauri star, and exhibits high X-ray flux, 
large lithium abundance, and evidence for an actively accreting disk with accretion rate 
estimated to be from 10~^ to 10~^ Mq yr~^ based on the large Ra equivalent width and 
the X-ray emission (Kastner et al. 2002 and refs therein). Despite its apparently advanced 
age, estimated to be 5 to 20 Myr, TW Hya is surrounded by a disk of mass ^ 5 x 10~^ Mq 
inferred from dust emission (Wilner et al. 2000). The disk is viewed nearly face-on, and is 
visible in scattered light to a radius of at least 3^.'5, or 200 AU at a distance of 56 pc (Krist 
et al. 2000; Trilling et al. 2001; Weinberger et al. 2002; Calvet et al. 2002). 

Several molecular species (CO, HCN, CN, HCO+, DCO+) have been detected in the TW 
Hya disk at millimeter and submillimeter wavelengths using single dish telescopes (Kastner 
et al. 1997; van Zadelhoff et al. 2001; Thi et al. 2001; van Dishoeck et al. 2003). The proximity, 
isolation, and rich chemistry of the TW Hya disk make this system an excellent target for the 
study of the physical and chemical structure of protoplanetary environments at high spatial 
resolution using interferometry. Models may be directly compared to aperture synthesis 
images of molecular lines to verify the inferences from low resolution observations, and to 
further probe the kinematics, physical conditions, and chemistry. A first step in this direction 
was the imaging analysis of HCO^ J=l-0 emission observed with the Australia Telescope 
Compact Array by Wilner et al. (2003), though the spatial and spectral resolution were 
insufficient to examine the kinematics in any detail. 

The prospects are excellent for interferometric observations of protoplanetary disks at 
submillimeter wavelengths (Blake 2003). The Submillimeter Array (SMA)-^ , now nearing 
completion on Mauna Kea, is an ideal instrument for making these observations (Ho et al. 
2004). In this Letter^ we present early SMA observations of the TW Hya disk in the CO 
J=3-2 and J=2^1 lines and dust continuum, with a maximum resolution of 2'^ The resulting 
images unambiguously reveal the rotation of the disk around the star and constrain the disk 
inclination and size. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2. 



Observations 



The SMA observations of TW Hya were made between 2003 March and June using two 
configurations of five of the 6 meter diameter antennas at each of two frequencies (230 GHz 
and 345 GHz). These observations provided 17 independent basehnes ranging in length from 
16 to 120 meters. Table 1 summarizes the observational parameters. The synthesized beam 
sizes were ^2'' at 345 GHz and ^4'' at 230 GHz with robust weighting. The SMA digital 
correlator was configured with both a narrow band of 512 channels over 104 MHz, which 
provided 0.2 MHz frequency resolution, or 0.18 km s~^ velocity resolution at 345 GHz, and 
several broader bands that provided 656 MHz of bandwidth for continuum measurements. 
Calibration of the visibility phases and amplitudes was achieved with observations of the 
quasar 3C273, typically at intervals of 30 minutes. Two weaker quasars closer to TW Hya 
on the sky, 1058+015 and 1037-295, were also observed to verify the effectiveness of phase 
referencing between 3C273 and TW Hya. Observations of Callisto provided the absolute 
scale for the flux density calibration. The uncertainties in the flux scale are estimated to 
be 20% based on the scatter of the measurements on different days, which is due to the 
uncertainties in pointing models and system temperature measurements. The data were 
calibrated using the MIR software package adapted for the SMA from the IDL version of the 
MMA software package developed originally for OVRO (Scoville et al. 1993). The continuum 
emission from TW Hya is strong enough for self-calibration, and an iteration of phase-only 
self-calibration was performed to improve the images by 20%. The NRAO AIPS package 
was used to generate continuum and spectral line images, all of which were CLEANed in the 
usual way. 



Calvet et al. (2002) have developed a physically self-consistent irradiated accretion disk 
model for TW Hya. This unique model requires a distinct transition in disk properties at 
a radius of ?^4 AU, likely due to a developing gap, and grains that have grown to sizes 
of ~1-10 cm in order to match the full spectral energy distribution from optical to radio 
wavelengths. The Calvet model is consistent with the radial brightness distribution of the 
disk observed at 43 GHz, where the dust emission has been resolved by the VLA (Wilner 
et al. 2000). The new SMA continuum measurements agree well with previous single dish 
observations at these wavelengths, and the predictions of the Calvet model. To compare the 
dust emission to the model radial brightness, we show in Figure 1 the visibility amplitudes at 
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Results 



3.1. 



Continuum Emission 
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345 GHz, annularly averaged in 20 kA bins, with the model curve (scaled by 0.9, well within 
the uncertainties of the flux scale). The good match of the SMA continuum measurements 
with previous data and the irradiated accretion disk model predictions gives confldence in 
the reliability of the submillimeter measurements, as well as the accuracy of the absolute 
flux scale. 



3.2. CO J=2-l and J=3-2 Emission 

Figure 2 shows spectra of the CO J=2-l and J=3-2 hne emission at the TW Hya 
continuum peak. Both lines are detected with high signal-to-noise ratios. The CO emission 
from the disk is spatially resolved; a circular Gaussian fit to the CO 3-2 integrated intensity 
map gives a FWHM size of 3'.'5 ± 0'.'3. Table 1 lists the integrated intensities and results 
of Gaussian fits to the central spectra. As indicated by Beckwith & Sargent (1993) and 
van Zadelhoff et al. (2001), CO J=2-l and CO J=3-2 emission are expected to be highly 
optically thick in circumstellar disks. These lines are thus excellent probes of the outer disk 
velocity field but do not robustly trace the disk mass (Simon et al. 2001; Qi et al. 2003). 
Further, the peak intensities should be tied to the kinetic temperatures of the regions of the 
disk where the lines become optically thick, and as TW Hya appears almost face on, any 
differences in brightness may be related to vertical temperature gradients. 

Comparing the intensities of the two transitions within the same region (here the CO 
J=3 2 data set is convolved with the CO J=2-l beam), the peak brightness of the CO 
J=3 2 line is 21.0 ± 0.7 K and that of the CO J=2-l line is 14.7 ± 0.4 K. Figure 3 shows 
the location of the r = 1 surfaces for these CO transitions calculated together with contours 
from the kinetic temperature model of Calvet et al. (2002) for three different depletion 
scenarios. All assume a constant depletion factor (Dc) of 10. The top set is calculated with 
the constant the one below includes additional jump depletion (Dj) by a factor of 10 
where the kinetic temperature falls below 22 K (the critical temperature of CO desorption, 
van Zadelhoff et al. 2001), and the bottom one includes of 100 at that same temperature 
threshold. For the first case, the temperature gradient probed by the two CO lines is very 
small, as both transitions are formed at nearly the same location in the disk. In the case 
where Dj is 10, however, the r = 1 surfaces for both transitions show an abrupt change 
near 130-150 AU in radius and 20-40 AU in height. With an even larger Dj of 100, the CO 
J=3-2 line becomes optically thick significantly deeper (colder) in the disk than does the 
CO J=2-l line in the region beyond 200 AU where the temperature is very low. The T^ of 
CO 3-2 is predicted to be smaller than that of CO 2-1 in this case, which is not compatible 
with the SMA spectra reported here. 
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The observed difference in the brightness temperatures of the CO hnes in Figure 2 must 
therefore arise from a combination of factors. The beam size for the CO J=3-2 spectrum 
corresponds to -160x120 AU (FWHM), that for CO J=2-l is -260x150 AU (FWHM). 
Both zones have a r = 1 surface above 22 K, though the 2-1 data do sample colder regions 
where the r = 1 surface drops to <20 K. The velocity-integrated intensities for the full disk 
in the two CO lines are 37.2 ± 7.5 Jy km s-^ for CO 3-2 and 12.4 ± 1.0 Jy km s-^ for CO 
2-1, as compared to previous single dish observations of 36.6 Jy km s~^ for CO 3-2 (van 
Zadelhoff et al. 2001) and 17.7 Jy km s'^ for CO 2-1 (Kastner et al. 1997). We regard these 
as consistent given the uncertainties in the flux scales of the various telescopes. The ratio 
of integrated intensities from the aggregate 345/230 GHz data sets is 2.45 ± 0.6, consistent 
with the value of 2.25 expected if the lines were optically thick and tracing isothermal gas. 
As mentioned above, even with the same convolving beam the CO J=3-2 emission remains 
brighter than that from CO J=2-l (21 K vs 15 K). Figure 2 shows that there is likely to be 
a small vertical temperature gradient to which these two transitions can be quite sensitive, 
e.g. for the case of an inclination angle of 7° where CO J=2-l traces gas near —15 K, while 
the more optically thick CO J=3-2 emission arises from gas closer to the disk surface and 
therefore somewhat warmer (—18 K). Given the likely flux uncertainties of the SMA and 
single dish observations, it is diflicult at present to quantify whether vertical temperature 
gradients in the disk or beam dilution dominate the observed difference in the brightness 
temperatures of the CO lines shown in Figure 2. Studies of additional CO transitions and, 
especially, isotopically substituted variants could settle this issue definitively. 

Figure 4 (upper panels) shows a set of velocity channel maps for 0.18 km s~^ intervals for 
CO J=3-2 hue. The main feature of interest is the position shift along a northwest-southeast 
axis at velocities that bracket the systemic velocity. In order to interpret these images in a 
quantitative way, we use a two-dimensional accelerated Monte Carlo model (Hogerheijde & 
van der Tak 2000) to calculate the radiative transfer and molecular excitation. We adopt the 
physical density and temperature structure derived by Calvet et al. (2002), which does an 
excellent job reproducing the dust continuum emission, and produce a grid of models with 
a range of inclinations and various depletion scheme to simulate the disk as imaged by a 
telescope with the resolution constrained only by the grid sampling (typically this is of order 
5-10 AU in the outer disk, or O'.'l — 0'.'2). Simulated observations of the model disks were 
produced by the MIRIAD software package using the UVMODEL routine to select synthetic 
visibility observations at the observed (^i, v) spacings. These synthetic visibilities were then 
processed into images in a manner identical to the SMA data for detailed comparisons (more 
detailed explanations and a description of the model may be found in Kessler 2003, for an 
application to LkCa 15 see Qi et al. 2003). 

The CO data provide a strong constraint on the disk inclination z, or more specifically sin 
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i. Previous analyses of scattered light images (Krist et al. 2000; Trilling et al. 2001) indicated 
that TW Hya and its disk are seen close to pole-on, but the inclination angle could not be 
well constrained. For the disk models used here, we fix the inner radius at 4 AU and the 
outer radius at 196 AU, and adopt Keplerian rotation along with a low turbulent velocity of 
0.05 km s~^ through comparison to the observed CO J=3-2 and J=2-l line profiles. Indeed 
turbulent velocities in the outer disk also affect the line width, but values of >0.1 km s~-^ 
or larger cannot reproduce the peanut-shaped morphology of the systemic velocity channel 
map while the line widths are essentially unchanged for turbulent velocities in the 0.05-0.1 
km s~^ range. Figure 2 shows the predictions for a range of inclination angles, which strongly 
impact the widths of the lines; the best fit inchnation is 7°, which is slightly larger than the 
inclination angle of 4° derived from the ellipticity of the infrared image of the disk around TW 
Hya in scattered light (Weinberger et al. 2002) where a circular disk is assumed. Specifically, 
the predicted line widths change by 30%, from 0.72 to 1.06 km s~^ for CO J=3-2 and from 
0.6-0.85 km s~^ for CO J=2-l, as the inclination angles vary from 6 to 8 degrees. Such 
significant line width variations with inclination angle (by as little as 1 degree) are expected 
when the disk is so close to pole-on because the locii of constant projected velocities is quite 
sensitivite to sin i - the error of which translates to a very tight limit on i. An inclination 
angle of 7° ± 1° is thus a robust fit that is strongly constrained by the CO line widths within 
the confines of the disk model employed. In general, the derived inclination angle can also 
be coupled to the size of the disk and the radial and vertical temperature gradients. The 
line intensity of CO 3-2 is still much larger than the model prediction, suggesting that the 
near-surface gas temperature in the Calvet model is somewhat too low. The line profiles are 
also not perfectly characterized as Gaussian (though they are very narrow), which we rely 
on for determining the line widths. Each these systematic uncertainties could be larger than 
the uncertainty of the inclination angle (1 degree) that we determine from the line widths, 
but can be constrained in the future by higher spatial resolution observations of additional 
CO transitions and isotopologues. 

A comparison of the observed channel maps of CO J=3-2 and synthetic maps generated 
by a detailed model with parameters listed in Table 2 is shown in Figure 4. While the model 
reproduces all of the major features observed, in particular the velocity gradient and peanut- 
shaped morphology of the central channel, the predicted CO 3-2 intensities in Figures 2 and 
4 are 40% smaller than observed. It is likely that vertical gas and dust temperature gradient 
difference is not quite captured by the SED-constrained disk models, which we will pursue 
further in a future paper. Resolved observations of CO isotopologues toward another T Tauri 
star, DM Tau, have been presented by Dartois et al. (2003), who used IRAM inteferometric 
data to constrain the depletion and temperature versus radius, as the less abundant species 
probe the inner part of the disk, within an LTE excitation/radiative transfer model. Such 
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an approach is useful for molecules such as CO that possess low critical densities, but as has 
been discussed by van Zadelhoff et al. (2003), LTE excitation is not a good approximation for 
higher dipole moment species such as CN. With the 2D Monte Carlo model, the statistical 
nature of the molecular excitation can be investigated within the context of realistic models 
of the disk structure. Chemical differentiation can therefore be constrained by resolved 
observations of multiple transitions of species with varying characteristics, and our future 
publications will thus include studies of other molecules, including CN and HCN, toward 
TW Hya. 

We thank all the SMA staff members for their diligent work in completing the SMA. 
We also thank P. D'Alessio and N. Calvet for providing the disk model of TW Hya and J.E. 
Kessler-Silacci for helpful discussions concerning the radiative transfer model. And we are 
grateful to J. H. Kastner for useful comments. 
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Table 1. Observational Parameters for SMA TW Hya 







ro 9-1 


Rest Frequency: 


345.796 GHz 


230.538 GHz 


Observations: 


2003 Mar 15 


2003 Mar 24 




2003 Apr 21 


2003 Apr 20 






2003 May 03 


Synthesized beam: 


2'.'9 X 2'/l PA 3.4° 


4'.'7 X 2'.'7 PA 9.9° 


R.M.S^ (continuum): 


35 mJy/beam 


6.8 mJy/beam 


Dust flux: 


1.46 ± 0.04Jy 


0.57 ± 0.02Jy 


Channel spacing: 


0.18 km s"^ 


0.26 km s"^ 


R.M.S.^ (hue): 


1.1 Jy/beam 


0.28 Jy/beam 


Integrated intensity 


23.7 K km s'^ 


10.6 K km s-i 


Peak intensity 


25.9 K 


14.7 K 




0.86 ±0.02 km s'^ 


0.73 ± 0.02 km s"^ 



^SNR hmited by the dynamic range. 

^Based on Gaussian fits to the spectra (FWHM). 



Table 2. Model Parameters Used in Simulating TW Hya CO emission 



Parameters 



Physical Structure Irradiated accretion disk (Calvet et al. 2002) 

Stellar Mass 0.6 M© 

Disk Size R^^ 4 AU, V^out^edge 196 AU 

Disk PA 45° 

Inclination Angle 7° 

Turbulent Velocity 0.05 km s~^ 

Depletion Factor 10 x (100 x for T<22 K) 
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Fig. 1. — Visibility amplitudes of the observed and modeled 345 GHz continuum emission 
from TW Hya, annularly averaged in 20 kA bins. The error bars represent ±1 standard 
deviation for each bin. The curve shows the visibility amplitudes derived from the model 
calculations of Calvet et al. (2002) of an irradiated accretion disk, which fits the full spectral 
energy distribution well after applying an amplitude scale factor of 0.9. 
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Fig. 2. — The CO J=3 2 and 3—2-1 spectra at the continuum (stellar) position, beam sizes 
are listed in Table 1. The curves show synthetic spectra for a small range of inclination 
angles using the model parameters listed in Table 2. 
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Fig. 3. — The solid contours show the gas kinetic temperature from the TW Hya model of 
Calvet et al. (2002). The dashed and dot-dashed lines show the r = 1 surfaces for the CO 
3—3-2 and 3—2-1 lines for the three depletion scenarios described in the text. 
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Fig. 4. — SMA images of the CO J=3-2 emission from TW Hya in 0.18 km s~^ velocity bins 
observed with the SMA {upper panels^ with a contour spacing of 5 Jy beam~^), together with 
images from a simulated observation using the model parameters listed in Table 2 (lower 
panels^ with a contour spacing of 3 Jy beam~^). The small cross indicates the position of 
the continuum source. 



